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ABSTRACT 


HexCycleSpanner (HCS) is a novel tool designed to achieve an elementary cycle refinement 
operation (ECRO) that is used to tighten a Directed Hamiltonian Circuit (dHC) - to be used 
iteratively as an extremely powerful technique to transform any given dHC to an optimum 
Shortest Directed Hamiltonian Circuit. Application to solving Asymmetric Travelling Salesman 
Problem (aTSP) is quite evident. ECRO uses a HCS to achieve an exchange of an outgoing arc- 
triplet of the given dHC with a corresponding incoming arc-triplet from outside the given dHC, 
thus resulting in a reconstructed transformed permuted dHC, while retaining the original 
direction/orientation of the given dHC. Effectively, the removal of the outgoing arc-triplet 
separates the given dHC into three cut-segments that are again rejoined by the incoming arc- 
triplet resulting in the reconstructed transformed permuted dHC wherein the original orientation 
of the dHC is retained - as defined by the orientation of the three cut-segments. This ECRO is 
analogous to the simplex pivot operation of Linear Programming wherein an exchange of the 
outgoing basic variable row with the incoming nonbasic variable column is performed to 
transform a given simplex tableaux representation into another equivalent simplex tableau 
representation that may be closer to the optimum tableaux. 


The concept of HexCycleSpanner can be generalized and extended to include incoming directed 
paths rather than incoming arcs to give rise to what may be called a CycleExpander (CyclExp) 
that may be used as an effective tool to construct a Directed Hamiltonian Circuit. 


The absence of a ‘global effectiveness measure’ (‘gem’) for aTSP - unlike in the case of LP - 
requires one to use a reference data base of optimum/shortest directed path between pairs of 
nodes to enhance the computational efficiency of the proposed iterative ECRO algorithm. 
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1. INTRODUCTION 


The earliest avatars of Hamiltonian Path and Hamiltonian Circuit can be traced all the way back 
to the earliest times of human intellectual endeavors, referring to the classic idea of the Knight’s 
tours on a chess board. Shortest (minimum weight) Hamiltonian Circuit (SHC) solves the 
Travelling Salesman Problem (TSP) by their very definition. TSP is one of the most widely 
studied combinatorial optimization problems, for which any improved solution algorithm is 
considered to be a critical achievement in terms of the computational complexity associated with 
such problems. 


Hamiltonian paths and Hamiltonian circuits have been extensively studied and reported because 
of its wide spread applications in various practical problem domains. Although various reports 
can be found in literature about the necessary conditions for Hamiltonicity and similarly on the 
sufficient conditions, there has not been any report of a strong/tight necessary & sufficient 
condition thereof. A necessary condition for a connected graph G(V,E) with |V| = n> 3 vertices 
to be Hamiltonian is established by Dirac’s [D52] theorem which states that the minimum vertex 
degree be at least equal to n/2 that is half the number of vertices n. However, this turns out to be 
very weak necessary condition. Chvatal & Erdos [CE72] gives sufficient conditions for Hamilton 
Path and/or Hamiltonian Circuit, using the concepts of stability number (maximum stable set or 
maximum independent set), connectivity and bridges; later generalized by Kouider [K94]. A 
theorem of Bondy & Chvatal [BC76] gives the necessary & sufficient condition for Hamiltonicity 
of a simple graph, using the concept of the closure of a graph, although its utility seems to be 
somewhat limited except for theoretical analyses. Lin & Kernighan [LK73] developed a heuristic 
approach known as the 2-opt move based on twin-arc-exchange as was originally proposed by 
Flood [F56] to improve a given Hamiltonian Circuit. Gould [G91] [G03] [G14] gives an 
excellent survey of the developments in terms of the Hamiltonian problem on undirected graphs. 


In the case of directed graphs, Bondy & Murty [BM08] mentions that Redei [R34] showed, by 
induction arguments on the number of vertices, that every tournament has a directed Hamiltonian 
Path. But every tournament need not be Hamiltonian, as for example a transitive tournament. 
Camio (1959) proved that every nontrivial strong tournament has a directed Hamiltonian Circuit. 
Chen & Manalastas [CM83] showed that every finite strongly connected digraph of stability 2 has 
a Hamiltonian path; for which Bondy [B95] gives a short and elegant proof. Using the concept of 
directed ear decomposition Knuth [K74] shows that every strong digraph admits a coherent 
feedback arc set. Bessy & Thomasse [BT03] [BT04] showed that every stong digraph admits a 
coherent cyclic order. Sebo [S07], Iwata & Matsuda [IM07] report results related to the existence 
of coherent cyclic order in strong digraphs. 


However there has neither been any report regarding any tight necessary & sufficient conditions 
for the existence of Hamiltonicity nor has there been any computationally efficient algorithms 
reported for determining Hamiltonicity in case of general networks (directed graphs/multigraphs). 
Zhang [ZOO] [CJMZO00] mentions that Kannellakis-Papadimitriou [KP80] arc-exchange update 
for asymmetric TSP adapted from Lin-Kernighan [LK73] edge exchange strategy for symmetric 
TSP seems to be a promising heuristic approach for dHC. Ban-Jensen & Gutin [BG09] [BG18] 
gives an excellent exposition of the developments in the algorithms for directed graphs. Gutin & 
Punnen [GP04] presents the details on the various developments in solving aTSP. 


Here we will see how an elementary cycle refinement operation (ECRO) is developed as a novel 


method to tighten a given Directed Hamiltonian Circuit (dHC) using a HexCycleSpanner (HCS) 
designed as a tool for the purpose - to be used iteratively as a computationally efficient algorithm 
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to determine the shortest directed Hamiltonian circuit and thus solve the asymmetric travelling 
salesman problem (aTSP). Also explained herein is the novel concept of a CycleExpander 
(CyclExp) that can be applied in the construction of a dHC. 


2. HEX-CYCLE-SPANNER 


A HexCycleSpanner (HCS) designed as a tool for the purpose of defining an elementary cycle 
refinement operation (ECRO) that transforms a given Directed Hamiltonian Circuit (AHC) into 
another one with improved properties - for example with reduced weight / length - is a 6-cycle 
(HCS is not a directed circuit) with six arcs and having some special characteristic features : it 
forms a cycle of six directed arcs, three alternating arcs oriented in one way and belonging to the 
given dHC whereas the other three alternating arcs directed in the opposite sense and belonging 
outside of the given dHC. An ECRO uses a HCS to achieve an exchange of judiciously chosen 
outgoing arc-triplet with a corresponding incoming arc-triplet, resulting in a reconstructed 
transformed permuted dHC. The HCS effectively cuts the given dHC into three cut-segments by 
removing outgoing arc-triplet and rejoins the three cut-segments using the incoming arc-triplet - 
in such a way that the resulting reconstructed transformed permuted dHC has exactly the same 
orientation as the original orientation of the given dHC - the original cycle orientation being 
defined by the orientation of the three cut-segments of the given dHC. Note that one canNOT 
achieve such an orientation-retaining transformation by cutting the given dHC into two pieces - 
three cut-segments is the minimum requirement to retain the original orientation in the 
transformed dHC. Thus an ECRO achieves an exchange between the outgoing arc-triplet and the 
incoming arc-triplet - this arc-triplet-pair forms the required HexCycleSpanner - that is a 6-cycle 
wherein the outgoing & incoming arcs are in alternating positions and in opposing orientation 
with respect to one another. 


3. STANDARD TEMPLATE FOR A GENERIC HEX-CYCLE-SPANNER 


Effective application of ECRO requires an efficient method of finding an appropriate HCS from 
among the various possible choices. To facilitate the search for an appropriate HCS it is essential 
to develop a standard template for a generic HCS. For this purpose, consider a 9-cycle dHC 
(1,2,3,4,5,6,7,8,9) with outgoing arc-triplet (9,1) (3,4) & (6,7) and incoming arc-triplet (9,4) (6,1) 
& (3,7) so that the given dHC gets cut into three cut-segments (1,2,3) (4,5,6) & (7,8,9) by the 
removal of the three outgoing arcs and these three cut-segments get rejoined by the three 
incoming arcs, resulting in a reconstructed transformed permuted 9-cycle dHC (1,2,3,7,8,9,4,5,6) 
wherein the original orientation of the given 9-cycle dHC is retained - as defined by the 
orientation of the three cut-segments (1,2,3) (4,5,6) & (7,8,9). 


Here, the HCS is the 6-cycle (9,1,6,7,3,4) - wherein the three outgoing arcs and the three 
incoming arcs are positioned mutually adjacent to one another - each outgoing arc is flanked by 
two incoming arcs on either end and each incoming arc is flanked by two outgoing arc on either 
end - the incoming arc-triplet oriented in the opposite sense relative to the orientation of the 
outgoing arc-triplet. Each outgoing arc of a given dHC has a one-to-one correspondence with a 
directed path from its tail-end to its head-end, containing the corresponding cut-segment of the 
given dHC flanked on either end by the same two incoming arcs that are adjacent to that outgoing 
arc, thus defining the cut-segment associated with this outgoing arc - this is an essential 
characteristic of a HCS. 


PIPR:©: Dr (Prof) Keshava. Prasad.Halemane Page 3 of 10 CycleExpander&HexCycleSpanner4aTSP 


PIPR:©: Dr (Prof) Keshava. Prasad.Halemane Page 4 of 10 CycleExpanderéHexCycleSpanner4aTSP 


Note that each of the interior nodes of these three cut-segments - represented by the nodes 2, 5, 8 
of the given dHC - can in fact represent a directed path joining the end nodes of the 
corresponding cut-segment, so that this template becomes a generic template applicable to any 
given dHC in general. 


As said earlier, an outgoing (og9/Arc) arc o9/ is flanked by two incoming arcs, one is an 
incoming (tail94Link) arc i94 incident on the tail-end of og9/Arc o09/ and points to the starting 
node of a cut-segment cs456 whereas the other is an incoming (head6/Link) arc i6] incident on 
the head-end of og9/Arc 09/ and links to the terminal node of this cut-segment cs456 - forming a 
directed path i94cs456i6/ from the tail-end of og9/Arc 09] to the head-end of the 09/ passing 
through (1) tail94Link i94 (2) cut-segment cs456 (3) head61Link i61 - establishing a one-to-one 
correspondence between the outgoing arc 09/ and the cut-segment cs456 - this cut-segment itself 
being flanked by these two incoming arcs, one being the tail94Link i94 and the other being the 
head61 Link i6/ forming the directed path i94cs456i6/. Observe similar association in case of the 
second outgoing arc og34Arc 034 with the tail37Link i37 the cut-segment cs789 and the 
head94Link i94 to form the directed path i37cs789i94 and also in case of the third outgoing arc 
0g67Arc 067 with the tail61Link i6] the cut-segment cs123 and the head37Link i37 to form the 
directed path i6/cs123i37. It is essential to note that the role of an incoming arc being either a 
tail-Link or a head-Link is relative to the specific outgoing arc og-Arc and not a permanent / static 
role. Also note that each of the three outgoing arcs get flanked by the very same two incoming 
arcs that flank the corresponding cut-segment to form the directed path from its tail-end to its 
head-end - this is an essential characteristic of a HCS. 


4. FINDING A HEX-CYCLE-SPANNER FOR A dHC 


Figure-1 gives the node arc incidence matrix for a generic 9-cycle dHC. Find an outgoing arc 
j1=@9,i1) which is flanked on each end by two incoming arcs, j10=(19,i4) at the tail-end and 
j11=(6,11) at the head-end, both being oriented in the opposite sense to that of this outgoing arc, 
and the farther end 14 & 16 of these two incoming arcs define a cut-segment (14,15,16) of the given 
dHC that has a one-to-one correspondence with the outgoing arc (i9,i1); the ends of this cut- 
segment being adjacent to another two outgoing arcs, j4=(13,14) and j7=(16,i7) that are joined by 
an incoming arc j12=(i3,i7). This set of three outgoing arcs (19,11) (3,14) & (16,17) along with the 
corresponding three incoming arcs (19,14) (16,11) & (43,17) form a 6-cycle HCS (19,11,16,17,13,i4). 
Of course, for each possible choice of an outgoing arc, there may be several choices of incoming 
arcs and therefore several possible choices of HCS that may be found, and the one that best 
satisfies the requirements, say for example reducing the overall cost/weight/length of the dHC 
needs to be chosen for performing the desired elementary cycle refinement operation. This 
process of selecting the matching pairs of outgoing & incoming arc-triplets is analogous to the 
selection of a matching pair of pivot row (outgoing basic variable) & pivot column (incoming 
nonbasic variable) in simplex pivot selection for linear programming problems. 


Note that the evaluation of a potential candidate (19,11) as an outgoing arc includes the 
confirmation of the existence of corresponding adjacent arcs (19,14) & (16,11) which are 
themselves adjacent to two outgoing arcs (13,14) & (16,17) and also an incoming arc (13,17) such 
that the given dHC gets cut into three segments (i1,12,13) (14,15,16) & (17,18,19) by the removal of 
the outgoing arc-triplet (19,11) (43,14) & (16,17) and these three segments again get rejoined by the 
incoming arc-triplet (19,14) (6,11) & (i3,i17) thus resulting in the reconstructed transformed 
permuted dHC (i1,i2,13,17,18,19,14,15,16). In this case, the effective reduction in the overall 
cost/weight/length indicated by (6w1 + 3w7 + 9w4) < (Qwl + 3w4 + 6w7) justifies the ECRO 
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defined by this exchange of the outgoing arc-triplet (19,11) (13,14) & (16,17) with the incoming arc- 
triplet (46,11) (3,17) & (19,14) as shown in Figure-2. 


L | jl | j2 | j3 | §4 | J5 | j6 | j7 | j8 | j9 | j10 | jl | jl2 
il | -—| + = 


i9 | + —| + 


Figure-1: Node Arc incidence matrix with cycle orientation indicated by ordering of nodes 


i il i2 i3 i4 i5 16 i7 i8 i9 
il + lw2 
i2 + | 2w3 
i3 + | 3w4 3w7 
i4 + 4w5 
i5 + 5w6 
i6 | 6wl + | 6w7 
i7 + | 7w8 
i8 + | 8w9 
i9 | 9wl 9w4 + 


Figure-2: Node Node Adjacency Matrix with Weights 


It may be worth to note a word of caution here. In the process of finding a matching pair of an 
outgoing arc-triplet with an incoming arc-triplet, if one is not careful to keep a watch on the 
relative positioning of the three outgoing arcs and the corresponding three incoming arcs in 
relation to the orientation of the given dHC, it is easy to get into a situation where the three cut- 
segments when rejoined to yield a reconstructed transformed permuted dHC may possibly end up 
being three mutually arc-disjoint cycles rather than a single dHC. For example, in the above case 
of a given dHC with a 9-cycle, if the three outgoing arcs (19,11) (13,14) & (16,17) are matched with 
the correspond incoming arcs (19,17) (13,11) & (16,14) the result is three mutually arc-disjoint 3- 
cycles (11,12,13) (14,15,16) & (17,18,19) which is an unintended result. 


To overcome this problem, or rather to prevent such an unacceptable situation, it is recommended 
to emphasize the fact that each of the three outgoing arcs, flanked in-between two incoming arcs, 
has a one-to-one correspondence with a directed path from its tail-end to its head-end, wherein 
the first and the last arcs are the two incoming arcs flanking that outgoing arc and the thus flanked 
segment of the directed path is indeed the cut-segment of the dHC corresponding to this outgoing 
arc. The outgoing arc (19,i1) associated with the directed path (19,(14,15,i6),11) which holds the 
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segment (14,15,16) flanked by the very same two end arcs (19,i4) & (16,i1) that are flanking the 
outgoing arc (19,11). This is an essential characteristic of a HCS. 


Now we can see that the earlier situation with the outgoing arc (i9,i1) flanked by the incoming 
arcs (19,17) and (13,i1) defines a cut-segment (i3,14,15,16,17) which along with these two incoming 
arcs flanking this outgoing arc doesn’t result in a directed path from the tail-end i9 to the head- 
end il; whereas this cut-segment again gets cut by the removal of its two end arcs (13,14) & (6,17) 
as outgoing arcs thus shrinking the cut-segment to (i4,15,16) while having no linking arcs between 
this cut-segment and the other two cut-segments (i1,i2,i3) & (17,18,19) thus resulting in three 
mutually arc-disjoint cycles. 


5. CYCLIC PERMUTATION 


In addition to bringing about an exchange between the outgoing arc-triplet and the incoming arc- 
triplet, the above defined ECRO using HCS can in fact be considered as the basic / elementary 
swapping operation in the process of generating cyclic permutations of dHC, wherein each ECRO 
brings about a single swap/exchange of the two of the three selected cut-segments w.r.t. the other 
in the cyclic order - as in the earlier example, transforming the given dHC (1,2,3,4,5,6,7,8,9) to 
the reconstructed transformed permuted dHC (1,2,3,7,8,9,4,5,6) wherein the two cut-segments, 
say (4,5,6) and (7,8,9) get swapped w.r.t. the third cut-segment that is (1,2,3). 


It is advised to keep track of the information associated with each of the potential candidate 
outgoing arc, as to the arc-triplet pair defining the HCS involved in the exchange as well as the 
change in the cost/weight/length associated with the corresponding ECRO. This will facilitate the 
determination of the best choice of the HCS and hence the best ECRO for the given dHC under 
consideration. 


6. CYCLE-EXPANDER 


The concept of HexCycleSpanner can be generalized and extended to define a CycleExpander 
(CyclExp) by allowing incoming directed paths instead of incoming arcs in the definition of a 
HexCycleSpanner. This will enable a CyclExp to be used as an effective tool in expanding a 
given directed circuit to include more and more nodes so as to result in the construction of a 
Directed Hamiltonian Circuit. 


7. CYCLE-EXPANDER AND HEX-CYCLE-SPANNER TO SOLVE ATSP 


The above concepts of CycleExpander (CyclExp) and HexCycleSpanner (HCS) can directly be 
applied in the construction of a dHC and on which ECRO can be performed using HCS in 
achieving further improvements in the dHC to determine the optimum/shortest directed 
Hamiltonian circuit thus solving the Asymmetric Travelling Salesman Problem - even for general 
networks / directed multigraphs. 


Although ECRO using HCS is analogous to simplex pivoting operation of linear programming 
there is no ‘global effective measure’ (‘gem’) defined for aTSP - unlike in the case for linear 
programming for which the ‘spdspds’ [H22] algorithm exploits the ‘gem’ defined by the 
‘infeasibility index’ to design a computationally efficient algorithm. In the case of aTSP the best 
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that may be done if at all is to maintain a reference data base of ‘optimum/shortest paths’ 
between pairs of nodes, which can be used effectively to retain such optimum path segments in a 
dHC when choosing the best possible HCS for the application of ECRO. 


We recommend that the Painted Network Algorithm of Rockafellar [R98] be adapted for path 
finding in the application of CyclExp and HCS to develop a computationally efficient solution 
strategy for aTSP. Appropriate fine-tuning of the algorithm with the use of CyclExp and HCS for 
the purpose of solving aTSP can result in extraordinary gains in computational efficiency. 


We recommend that one initiates a directed ear decomposition [BG09] of the given network 
(directed graph / multigraph) using a depth first search on minimum cost outgoing arc at every 
node and to choose - either the primary ear circuit defined here as the one wherein every arc 
corresponds to the optimum (minimum cost/weight/length) outgoing arc from its initial node - or 
the leap ear circuit defined here as the one obtained by conjoining the primary ear circuit with 
other suitable ones through a segment exchange operation to form the maximal cardinality 
(number of nodes) ear while possibly compromising on the cost/weight/length associated with the 
arcs at the conjoining nodes - as a good starting point to apply the CyclExp for constructing a 
dHC which can then be further refined / improved through ECRO by the application of HCS. 


One may choose the best possible HCS that will potentially result in the best possible reduction in 
the cost / weight / length as a result of the application of the corresponding ECRO. For this 
purpose, it may be useful to keep track of the data associated with every potential choice of HCS 
defined by the corresponding pair of arc-triplets, and choose the best among them that result in 
the best possible improvement in the reconstructed / transformed / permuted dHC. 


Alternatively, one may consider some pre-processing using the concept of bead-bridge 
decomposition into strong components may provide useful information to possibly enhance the 
computational efficiency by working on the components and then integrating the results later. 


Also, some appropriate shortest path algorithm (e.g. Floyd-Warshall) may be applied on the given 
network, initially, in order to develop the required reference data base of min-cost directed path 
segments between node pairs, which will later facilitate in the decision regarding the best possible 
choice for a HCS. If any sub-sequence of any of these paths is a matching cycle-segment of the 
dHC under investigation, then such a cycle-segment cannot be subjected to any further cuts by 
application of ECRO using any HCS, because it corresponds to the best possible choice in terms 
of minimum-cost / least-weight / shortest-length sub-sequence / cycle-segment of the dHC and 
hence no further improvement can possibly be achieved through the application of any ECRO 
using whatever available HCS. 


Alternatively, one may even explore the possibility of starting with the optimum arborescence 
along with the corresponding root-location problem mentioned by Fujishige & Kamiyama [FK12] 
[K14] instead of the Floyd-Warshall type approach, before the actual application of the CyclExp 
to construct a dHC and then proceed with the application of ECRO with HCS. 
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8. CONCLUSION 


The ECRO operation and the associated tools CyclExp & HCS facilitate the design of algorithms 
for constructing dHC for a given network (directed graph / multigraph) and also to refine a given 
dHC to determine the optimum (minimum-cost / least-weight / shortest-length) dHC and hence 
solve aTSP. ECRO using CyclExp & HCS is flexible to adopt various levels of fine-tuning that 
may be required, irrespective of the details of whichever algorithm is used, in order to address 
various possible requirements including one of computational complexity. 
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